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ABSTRACT 

A new efficient scheme for imaging gamma-ray emitting ob¬ 
jects is advocated in this work. A few years back we have 
put forward the idea of collecting data on Compton scattered 
rays to reconstruct an object in three dimensions and shown 
how it works with a collimated detector. To increase drasti¬ 
cally the sensitivity of this modality as well as its field of view 
and resolution we propose now that data acquisition should be 
performed without collimation. We discuss image formation 
by scattered radiation in this context by computing and com¬ 
paring the related Point Spread Functions (PSF) and comment 
on their properties. We also present numerical simulations to 
support the attractiveness of this modality. 

Key words: Gamma-ray imaging, sensitivity, photon Comp¬ 
ton scattering, image formation, object reconstruction. 

1. INTRODUCTION 

Emission imaging with gamma rays is widely used in numer¬ 
ous fields such as medical imaging, non-destructive testing, 
gamma astronomy and environmental survey. In conventional 
nuclear imaging, a collimated gamma camera rotates in space 
to collect primary radiation emitted by an object under in¬ 
vestigation. In this case Compton scatter radiation behaves 
generally as noise hindering image quality and consequently 
correction to scatter should be applied. 

However recently an interesting new imaging concept, which 
precisely uses, as imaging agent, scattered radiation by the 
object medium (instead of primary radiation), has been pro¬ 
posed. A spatially fixed collimated gamma camera records 
now images labeled by the energy of scattered radiation (or 
equivalently its scattering angle). It is then shown that the 
reconstruction of a three-dimensional object is feasible us¬ 
ing this data [1, 2, 3], However in this situation, the im¬ 
age sensitivity is considerably affected due to the presence 
of the collimator. Only about one out of 10 4 scattered pho¬ 
tons reaches the detector. Therefore in order to record a much 
larger amount of scattered radiation, we propose to extend the 


working principle of this gamma-ray imaging to a functioning 
modality without collimator, as depicted by Fig. 1. However, 
as can be seen, this procedure differs from other proposals on 
scattered radiation imaging processes [4, 5, 6, 7], in particu¬ 
lar: 

- Compton tomography [8], which reconstructs the elec¬ 
tron density of the object (instead of its activity density), and 
uses a moving point-like detector collecting scattered radia¬ 
tion from an external radiation source, 

- Compton camera [9], which reconstructs the activity den¬ 
sity of an object from scattered radiation using coincidence 
measurements between a site on a scatter layer-detector and 
another site on an absorption-detector. 

As the true three-dimensional problem implies involved 
space geometry of scattered rays, we shall first study its two- 
dimensional counterpart to test the viability of this idea. To 
this end, we give a careful analysis of the image formation. 
The corresponding PSF is derived and shall be compared with 
the previous case where a collimator is present. Next we per¬ 
form numerical reconstruction of the object radiation activ¬ 
ity from simulated scattered radiation (recorded by a camera 
without collimator) to illustrate the working of this concept. 



Detector 


Fig. l. Two modalities in emission scattered radiation imag¬ 
ing. 






2. IMAGE FORMATION 


To understand image formation by scattered radiation, we fol¬ 
low radiation propagation in a two-dimensional (2D) slice (or 
thin section) of a scattering medium with an electronic den¬ 
sity n e assumed to be approximately constant. This can be 
in principle implemented with a slit placed above a detec¬ 
tor without collimator, see Fig. 2. Moreover, to concentrate 
chiefly on scattering effects radiation, attenuation will be left 
out as working hypothesis. 



Fig. 2. Two dimensional working with a slit. 


Let 

- S = (xs, ys ) be a radiation emitting point source of an ob¬ 
ject described by its activity density /(S) (number of photons 
emitted per unit time isotropically in all directions, 

- M = (xm,Vm), a scattering site in the medium, 

- n e ( M), the electron density at site M, 

- D = (xd, Vd), a detection site on a linear detector which 
collects the photon flux density at a photon energy E u (see 
Eq. (1)). 

The energy of scattered photon is related to the scattering 
angle lu by the Compton relation: 


E u — Eq 
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( 1 ) 


where Eq is the emitted photon energy, e = Eq/ me 2 and me 2 
the rest energy of the electron. 


From Fig. 3, one can see that the photon flux density 
reaching a site M is the number of photons emitted into the 
angular fan dtt per unit length and per unit time: 
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where \SM\ is the distance between sites S and M and dS$ 
the area element around S. The fraction of photons scattered 



Fig. 3. Geometry of Compton scattering: source site S, scat¬ 
tering site M and detection site D of width a'. 


in the direction making an angle u> with the incident direc¬ 
tion depends on the Compton differential cross section <jg (w) 
(which has the dimension of a length in 2D) and on the num¬ 
ber of electrons n e (M) c?Sm at site M, dSw being the inte¬ 
gration area element around M. Hence the scattering photon 
flux density received at the detection site D is given by : 
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where 6 is the angle between the outgoing photon unit vector 
with the detector normal unit vector, \MD\ the distance from 
scattering site M to detection site D. In fact, for a given point 
source S, there will be two scattering sites Mi et M 2 located 
on two arcs of circle subtending a scattering angle (tt — uj), as 
shown in Fig. 4. 

The total photon flux density at a site D is g(D, u>), the 
integral over all source sites and all scattering sites such that 
the scattering angle is to. This last constraint is expressed by 
a (5-function in the integration as 


3(D,w) 


#(D,w|S) S(SMD- (tt-w)), (4) 


where SMD is the angle at vertex M of the triangle SMD 
(see Fig. 4). Eq. (4) describes the basic image formation pro¬ 
cess from scattered radiation measured on a detector without 
collimator. 


3. COMPUTATION OF THE IMAGE OF A POINT 
SOURCE 

We now compute and study the essential object in this new 
imaging concept: the Point Spread Function (PSF). It is, by 
definition, the image of a single point source at site S. Eq. (4) 
can be rewritten in terms of the PSF as follows : 
































9{ D,W)= J dSs /(S) PSF( D,w]S). (5) 

As shown before for a camera without collimator, the scat¬ 
tering sites due to a single point source S are located on two 
circular arcs subtending an angle ( 7 r — ce) (Eq. (4)). The 
photon flux density received at site D is then given by an in¬ 
tegration over these two circular arcs (see Fig.4). 



Fig. 4. Locus of scattering sites due to single point source 


For computation ease, we shall assume a uniform electron 
density n e and use polar coordinates (r, 7 ), such that S = 
(d, a) with DS = d and M = (r, 7 ), with DM = r and DM ■ 
DS = rd cos 7 . 

The circular arcs have polar equations: 
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f = <i sm(a, + 7 (6) 
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where 7 is the angle between DS and DM. The distance 
\SM\ can be extracted from a simple identity in the triangle 
DSM: 


\SM\ = d^ 


(7) 


sinw 

And the integration area <7 Sm is now reduced to the arc ele¬ 
ment: 


\J dr 2 + r 2 d"f 2 = —r— dj. ( 8 ) 

smw 

Hence the PSF is given by the sum of the two integrations 
on 7 : 
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where K(uj) = An e ag (ce)/o/27r, <jg{ui) the differential Comp¬ 
ton cross section at scattering angle w, fo the intensity of the 
single point source and cos 9 = sin(a — 7 ) if the detector 
lies along the Ox axis and l is the distance between the line 
detector and the linear lower boundary of the medium L. The 
integration is carried out over the points inside the scatter¬ 
ing medium. Therefore when the medium is of finite extent, 
the limit of the integration 7 i(u>), which corresponds to the 
intersection of the arcs of circle with the scattering medium, 
should be calculated beforehand, see Fig. 4. 



Fig. 5. Scattering site in the presence of a collimator from 
single point source 


Now if the collimator is mounted on the detector, then 
only one scattering site M, located on the perpendicular to 
the detector at site D, will contribute to detection site D, (see 
Fig. 5). Thus the integration on 7 is restricted by a delta 
function which picks out only the corresponding value of 7 , 
i.e.: 

7r 

Icol = 2 ~ a ' ( 10 ) 

The resulting PSF expression for a collimated detector is: 


PSF col ( D,w|S) 
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Now at a fixed scattering angle ce, the PSF curve as func¬ 
tion of the detector position with collimator has a Mexican 
hat shape (see Fig. 6 ) 






























Fig. 6. PSF shapes with collimator 


while the PSF curve without collimator has a wide Lorentzian 
shape (see Fig 7). 



Fig. 7. PSF shapes without collimator 

For example at a scattering angle of 30 degrees. Fig. 8 
shows that the PSF without collimator is about 10 times stronger 
than the PSF with collimator. 


4. NUMERICAL RECONSTRUCTION RESULTS 

As an illustration of this new imaging concept, we carried out 
numerical reconstructions of a two-dimensional Shepp-Logan 
medical phantom from data computed with our model. 

The 2D original object (see Fig. 9) is placed at the center 
of the scattering medium and a unit distance above the detec¬ 
tor. A line detector of 55 pixels of 1 unit length, placed on 


Scattering angle = 30 degrees 



Fig. 8 . Comparison of PSF with (lower blue line) and without 
collimator (upper red line) 

the axis y = 0, is simulated. We consider that the scatter¬ 
ing medium has the same properties as water. It consists of a 
discretized square of 55 x 55 elements of unit area. 



Fig. 9. Original Shepp-Logan phantom 


A series of 55 images of the object corresponding to 55 
different scattering angles (12° < u < 132°) have been sim¬ 
ulated. We construct the 3025 x 3025 weight matrix by com¬ 
puting, for each mesh point source, the PSF at the different 
scattering angles for each site on the detector. The recon¬ 
struction is carried out by inverting this weight matrix using 
the Singular Value Decomposition method, which is less time 



















5. CONCLUSION 


consuming, compared to other reconstruction methods. 


The feasibility of image reconstruction using Compton scat¬ 
tered rays detected by a gamma-camera without collimator, 
operating in a fixed position, is demonstrated in this study. 
This is the essence of a new concept of high sensitivity imag¬ 
ing, which takes advantage of scattering rays instead of re¬ 
jecting them as done usually. The main point in this imaging 
process by Compton scattered radiation is the fact that data 
acquisition is performed without the usual motion of the de¬ 
tector. This is a major advantage compared to existing imag¬ 
ing systems which require a heavy, bulky and costly mechani¬ 
cal rotation mechanism to move the detector around in space. 
Work towards an extension to three-dimensional imaging is in 
progress. The modeling and simulations of multiple scattered 
radiation in this context will be also subjects of future inves¬ 
tigations. These promising results may open the way to new 
high sensitivity imaging devices which will have applications 
in nuclear medicine, non-destructive industrial control, high 
Fig. 10. Shepp-Logan phantom reconstruction with collima- energy astrophysics, environmental survey, etc. 
tor 




Fig. 11. Shepp-Logan phantom reconstruction without colli¬ 
mator 

Fig. 10 and Fig. 11 show the reconstruction results with 
and without collimator. One observes a better agreement with 
the original object when the collimator is removed. In Fig. 
10 , we can see that the part of the object near the detector 
is better reconstructed than the upper part of the object. The 
three small structures are invisible. The reconstructive rel¬ 
ative error is about 9.13 %. But without collimator (see Fig. 
11), the whole object is correctly reconstructed. All structures 
are visible and the relative error is about 1.17 x 10 -4 %. 

These results have been subsequently validated by Monte- 
carlo simulations [10]. 
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